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Abstract 

We investigate the possible frictionless transport of many composite (condensed) fermions 
at room temperature regime along an annular tube with transversely wavy-corrugations by 
using the verified transition-rate model and boundary perturbation approach. We found 
that for certain activation volume and energy there exist possible frictionless states at room 
temperature regime. 
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1 Introduction 

Understanding the frictionless transport of system of (interacting or noninteracting) fermions is 
crucial to our knowledge of the nature considering the neutron star, the superconductivity, the 
superfluidity, etc. One relevant study is the existence and characteristics of the Fermi surface 
for a system of interacting fermions. Usually such a surface is introduced into momentum space 
only for a system of noninteracting fermions. 

To be specific, this surface represents the limit of occupation of the different single-particle mo- 
mentum states in the ground state of the system. All the states with momentum contained 
within this surface are occupied, all those with momentum outside this surface are unoccupied. 
However, Migdal [1] noticed that, under some circumstances, the mean occupation number of 
different single-particle momentum states in the true ground state still possessed a discontinuity 
for a system of interacting fermions. Meanwhile, as Luttinger pointed out, the existence of the 
Fermi surface depends on the nature of forces between the fermions [2]. 

Other relevant study is, e.g., thermodynamical properties of trapped noninteracting Fermi gases 
in gravitational fields [3] which was originated in view of the successful experiments with the 
trapping and cooling of weakly interacting fermionic isotopes [4]. There is no doubt that it is 
important to theoretically investigate the thermodynamical and statistical properties of system 
of fermions even there are many approaches [5-11]. Here we shall investigate the possible fric- 
tionless transport of many fermions in a confined domain via a different approach considering 
the shearing response of many-fermion system. Note that researchers have been interested in the 
question of how matter responds to an external mechanical load. External loads cause transport, 
in Newtonian or various types of non-Newtonian ways. Amorphous matter, composed of poly- 
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mers, metals, or ceramics, can deform under mechanical loads, and the nature of the response to 
loads often dictates the choice of matter in various applications. To the best knowledge of the 
author, the simplest model that makes a prediction for the rate and temperature dependence 
of shear yielding is the rate-state model of stress-biased thermal activation [12-14]. Structural 
rearrangement is associated with a single energy barrier E that is lowered or raised linearly by 
an applied stress a : R± = vq exp[— E/^ksT)] exp[ztaV* / (ksT)], where ks is the Boltzmann 
constant, vq is an attempt frequency and V* is a constant called the 'activation volume'. In 
amorphous matter, the transition rates are negligible at zero stress. Thus, at finite stress one 
needs to consider only the rate i?+ of transitions in the direction aided by stress. 
The linear dependence will always correctly describe small changes in the barrier height, since 
it is simply the first term in the Taylor expansion of the barrier height as a function of load. 
It is thus appropriate when the barrier height changes only slightly before the system escapes 
the local energy minimum. This situation occurs at higher temperatures; for example, New- 
tonian transport is obtained in the rate-state model in the limit where the system experiences 
only small changes in the barrier height before thermally escaping the energy minimum. As the 
temperature decreases, larger changes in the barrier height occur before the system escapes the 
energy minimum (giving rise to, for example, non-Newtonian transport). In this regime, the 
linear dependence is not necessarily appropriate, and can lead to inaccurate modeling. To be 
precise, at low shear rates (7 < j c ), the system behaves as a power law shear-thinning material 
while, at high shear rates, the stress varies affinely with the shear rate. These two regimes 
correspond to two stable branches of stationary states, for which data obtained by imposing 
either a or 7 exactly superpose. 

In this paper, considering the general forcing (either gravity [3] or electric field [6]) we shall 
adopt the verified transition-rate-state model [12-14] to study the frictionless transport of many 
fermions within a corrugated annular tube. To obtain the law of shear-thinning matter for ex- 
plaining the too rapid annealing at the earliest time, because the relaxation at the beginning 
was steeper than could be explained by the bimolecular law, a hyperbolic sine law between 
the shear (strain) rate : 7 and shear stress : r was proposed and the close agreement with 
experimental data was obtained. This model has sound physical foundation from the thermal 
activation process [12-14] (a kind of (quantum) tunneling which relates to the matter rearrang- 
ing by surmounting a potential energy barrier was discussed therein). With this model we can 
associate the (shear-thinning) fluid with the momentum transfer between neighboring atomic 
clusters on the microscopic scale and reveals the atomic interaction in the relaxation of flow with 
dissipation (the momentum transfer depends on the activation (shear) volume : V* = Vh which 
is associated with the center distance between atoms and is equal to fc^T/ro (T is temperature 
in Kelvin, and To a constant with the dimension of stress). 

To consider the more realistic but complicated boundary conditions in the walls of the annular 
tube, however, we will use the boundary perturbation technique [15-16] to handle the presumed 
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wavy-roughness along the walls of the annular tube. To obtain the analytical and approximate 
solutions, here, the roughness is only introduced in the radial or transverse direction. The rel- 
evant boundary conditions along the wavy-rough surfaces will be prescribed below. We shall 
describe our approach after this section : Introduction with the focus upon the transition-rate 
approach and boundary perturbation method. The approximate expression of the transport is 
then demonstrated at the end. Finally, we will illustrate our results and give discussions therein. 

2 Formulations 

We firstly take into account the condensed system of many fermions subjected to random ther- 
mal fluctuations (under external forcing). In thermally-activated motion, mobile fermions in 
a many-fermion system may interact with other fermions even they are already in a preferred 
motion. The rate of deformation (strain) is controlled by the rate at which thermal energy can 
help the composite systems overcome their energy barriers, allowing the rest of other fermions 
to spread. 

The attractive interactions involving pairs of moving fermions lower the energy as the mech- 
anisms that enable them to form a composite moving subsystem cost energy. Repulsive in- 
teractions, on the other hand, require less work to overcome, and do not usually transform 
the composite subsystem or leave residual subsystem after the interaction. Parts of condensed 
fermions easily pull away after having been forced to pass the repulsive obstacle by an exter- 
nal stress. Accordingly, attractive interactions require thermal activation and are said to be 
temperature-dependent, while repulsive ones are not. 

In fact, thermal energy is supplied by random thermal fluctuations, and motion of composite 
(condensed) fermions depends on the number of fluctuations that supply the interacting sub- 
systems the energy they need. The number of such successful outcomes is N s = A^P (success), 
where N is the number of attempts in the complete many-fermion system. The probability of 
success is the probability that the thermal jump Uj is greater than AE, the energy required 
to surmount the barrier and assumed to follow an Arrhenius law [12-14] which is P(success) = 
P(Uj > AE) = exp(— AE/ksT). Considering the strain gained at each successful attempt leads 
to the thermal- activation controlled expression for the strain rate 7 P = A f Vo exp(— AE{t s ) /UbT). 
Here 7 P0 is a intrinsic value that has units of strain rate and depends on the average wait- 
ing time at the intersection point, the strain released after the events, and the frequency of 
thermal fluctuations, which is some fraction of the Debye frequency (10 13 s -1 ). We remind 
the readers that AE(t s ) is a function of r s which is a concentrated shear stress due to the 
short-range interaction with another subsystem (or obstacle), emphasizing the localized nature 
of thermally-activated events. It is possible to attribute AE in above expression to one type of 
thermally- activated process. This simplification applies when one process has a much smaller 
AE than the rest. Thermally controlled deformation, however, is a complex collective phe- 



4 



nomenon of many thermally-activated processes. To continue to use this theory, AE must be 
treated as an effective energy covering all possible types of such processes. 
We shall consider a steady transport of many fermions in a wavy-rough annular tube of r\ 
(mean-averaged inner radius) with the inner interface being a fixed wavy-rough surface : r = 
r\ + esin(/c# + /3) and ri (mean-averaged outer radius) with the outer interface being a fixed 
wavy-rough surface : r = r 2 + esin(A;0), where e is the amplitude of the (wavy) roughness, /3 
is the phase shift between two walls, and the roughness wave number : k = 2ir/L (L is the 
wavelength of the surface modulation in transverse direction). 

Firstly, this matter (composed of many condensed (composite) fermions) can be expressed as 
[12-14] 7 = 70 sinh(r/ro), where 7 is the shear rate, r is the shear stress, tq = 2kBT/Vh, and 
7o(= CkksT exp(—AE / ksT) / h) is with the dimension of the shear rate; here Ck = 2Vh/V m is 
a constant relating rate of strain to the jump frequency (V^ = A2A3A, V m = A2A3A1, A2A3 is 
the cross-section of the transport unit on which the shear stress acts, A is the distance jumped 
on each relaxation, Ai is the perpendicular distance between two neighboring layers of particles 
sliding past each other), accounting for the interchain co-operation required, h is the Planck 
constant, AE is the activation energy. 

In fact, the force balance gives the shear stress at a radius r as r = —(rSQ)/2 [15]. SQ is the 
net effective external (gravity or electric field) forcing along the transport (or tube-axis : z-axis) 
direction (considering dz element). Introducing the forcing parameter (f> = — {t2/2tq)5Q then we 
have 7 = 7osinh(</>r/r2). As 7 = —du/dr (u is the velocity of the transport in the longitudinal 
(z-)direction of the annular (cosmic) string), after integration, we obtain 

u = u s + ^^[coshi^ — cosh(— )], (1) 

9 r 2 

here, u s (= u s u p ) is the velocity over the (inner or outer) surface of the annular (cosmic) string, 
which is determined by the boundary condition. We noticed that a general boundary condition 
for transport over an interface [15] was 

Su = L° 7 (l - ^)" 1/2 , (2) 

7c 

where 5u is the velocity jump over the interface, L° s is a constant slip length, j c is the critical shear 
rate at which the slip length diverges. The slip (velocity) boundary condition above (related to 
the slip length) is closely linked to the mean free path of the particles together with a geometry- 
dependent factor (it is the quantum-mechanical scattering of Bogoliubov quasiparticles which 
is responsible for the loss of transverse momentum transfer to the confined surfaces [17]). The 
value of 7 C is a function of the corrugation of interfacial energy. 

With the slip boundary condition [15], we can derive the velocity fields and transport rates 
along the wavy-rough annular tube below using the verified boundary perturbation technique 
[15-16] and dimensionless analysis. We firstly select L° s to be the characteristic length scale and 
set r' = r/L°, R\ = 77 /L°, R2 = ^/-L^, e' = e/L° s . After this, for simplicity, we drop all the 
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primes. It means, now, r, R±, R 2 and e become dimensionless {<p and 7 also follow). The wavy 
boundaries are prescribed as r = R 2 + e sm.(k9) and r = R\ + e sin(/c# + /3) and the presumed 
steady transport is along the z-direction (annulus-axis direction). 

2.1 Boundary Perturbation 

Along the outer boundary (the same treatment below could also be applied to the inner bound- 
ary), we have 7 = (du) / (dn)\ on interfaces- Here, n means the normal. Let u be expanded in e 

u = u + eui + e 2 u 2 H , 

and on the boundary, we expand u(ro + edr, 9(= 9q)) into 

dv 2 

u ( r ,9)\(r +edr,9 ) = u( r 0,&) + 4 dr u r (rQ, #)] + <? [ — U„ (r , 9)} H = 

{uslip + — [cosh0 - cosh( — )]}| on interfaces, r = R 1} R 2 ; (3) 

where 

u slip\on interfaces ~~ on interfaces' 

(4) 

7c 

Now, on the outer interface (cf. [16]) 

du _ V(r - R 2 — esm(k9)) M 2 ^ 2 o,. ml _i r 
7 = ^ = VU ■ \V(r-R 2 -esin(k9))\ = [1 + 6 7* C ° S ^ 2 KI(*+«M>- 

e ^2 cos (^WI(i?2+«fr,0)] = u 0r \R 2 +e[ui r \ R2 +u 0rr \n 2 sin(k6)- 



^u%\r 2 cos(k9)] + e 2 [ 

1 . 9„ 



— «oJfl 2 cos(fc0)] + e 2 [--— cos 2 (k9)u 0r \R 2 + u 2r \n 2 +u lrr \ R2 sm(k9)+ 



o u 0rrr \R 2 sin 2 (k9) - — cos(k9)(u lf) \r 2 + u 0f)r \r 2 sm(k9))\ + 0(e 3 ). (5) 

Considering L° ~ Ri,R 2 S> e case, we also presume sinh^ <C 7 c /7o- With equations (1) and 
(5), using the definition of 7, we can derive the velocity field (u) up to the second order : 

u (r,9) = -(R 2 j /cf)){cosh((f)r/R 2 ) - cosh^ [1 + e 2 <j) 2 sin 2 {k9) / {2R 2 2 )} + 

ecpsinhcp sin(k9)/R 2 } + u sHp \ r=R2+esinlyke) . 

The key point is to firstly obtain the slip velocity along the boundaries or surfaces. After lengthy 
mathematical manipulations, we obtain the velocity fields (up to the second order) and then we 
can integrate them with respect to the cross-section to get the transport (volume flow) rate (Q, 
also up to the second order here) : 

9 P rR2+£sin(ke) 
/0 /Ri+esin(fc0+/3) 



Q = / / u(r, 9)rdrd9 = Q + e Q Po + e 2 Q P2 . 

JO JRi+esm(k6+8) 
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In fact, the approximate (up to the second order) net transport (volume flow) rate reads : 

Q = 7T7o{L°(i?l - R\) sinh</> (1 - !^)-i/2 + ^[ (jR 2 _ tffi cosh0 _ h R \ sinh0 _ 

7c/7o <P <P 

i?ii? 2 sinh(<A) + ^(cosh^ - cosh(<A)]} + e 2 {^u s% > (#l " 
tii <P til & 

L^7osmh0(l + S ^)(-k 2 + «/, 2 )[l - A 2 ] + ^o^sinh^) - i? 2 sinh</>]- 

4 7c/ 70 -K2 ^ ^2 

7T i?2 r w -^1 \ l ^ ,,_ -R? , 

7770 — [cosh - cosh(</>— )] + -7o</>cosh</>[i? 2 - — ] + 

7T7o{[sinh(/) + L° s cosh 0(1 + ^^)](i? 2 - R\ cos/3)} + ^o^ 2 - cosh(/>+ 

7c/7o ^ </> 

r0 7T 2 . COsh</> /-Rl\2i 



L «i^v¥ [1 - ( ^ )llcosM ' (6) 



Here, 

^ o =^7o[sinh0(l-^)-^]. ( 7) 

7c/7o 

3 Results and Discussions 

With above detailed derivations, now, we firstly check the presumed wavy-roughness effect 
(or combination of curvature and confinement effects) upon the possible frictionless transport 
of many condensed (composite) fermions because there are no available experimental data and 
numerical simulations for the same geometric configuration (annular tube with wavy corrugations 
in transverse direction). With a series of forcings (due to externally imposed gravity or electric 
field forcings) : <fi = — i? 2 (5^/)/(2ro), we can determine the enhanced shear rates (drf/dt) due to 
these forcings. From equation (5), we have (up to the first order) 

= [sinh </> + e sm{k6) 4~ cosh <f>] . (8) 
at at I%2 

The parameters are fixed below (the orientation effect : s'm(k9) is fixed here). r 2 (the mean 

outer radius) is selected as the same as the slip length L^. The amplitude of wavy roughness 

can be tuned easily. The effect of wavy-roughness is significant once the forcing ((/>) is rather 

large (the maximum is of the order of magnitude of e[</>tanh(</>)/i? 2 ]). 

If we select a (fixed) temperature, then from the expression of To, we can obtain the shear stress 
r corresponding to above gravity forcings (<fi) : 

t = r sintr 1 [sinh(0) + e sin(fc0) cosh(</>)] . (9) 

There is no doubt that the orientation effect (9) is also present for the condensed many-fermion 
system. For illustration below, we only consider the maximum case : |sin(/c#)| = 1. The wave 
number of roughness in transverse direction is fixed to be a constant. 
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As the primary interest of present study is related to the possible frictionless transport or 
formation of superfluidity (presumed to be relevant to the many-fermion system as mentioned 
in Introduction) due to strong shearing, we shall present our main results in the following. Note 
that, based on the absolute-reaction-rate Eyring model (of stress-biased thermal activation), 
structural rearrangement is associated with a single energy barrier (height) AE that is lowered 
or raised linearly by a (shear) yield stress r. If the transition rate is proportional to the plastic 
(shear) strain rate (with a constant ratio : Co; 7 = CoRt, Rt is the transition rate in the direction 
aided by stress), we have 

where v is an attempt frequency or transition rate, CqVq ~ 70 exp(AE/ksT), or 

r = 2^f-^eM^E/k B T) if (11) 

V h CqUq k B T 

It is possible that the frictional resistance (or shear stress) can be almost zero (existence of 

r ~ 0) from above equations (say, equation (10) considering a sudden jump of the resistance). 

The nonlinear character only manifests itself when the magnitude of the applied stress times the 

activation volume becomes comparable or greater in magnitude than the thermal vibrational 

energy. 

Normally, the value of Vh is associated with a typical volume required for a microscopic shear 
rearrangement. Thus, the nonzero transport rate (of the condensed many-fermion system) as 
forcing is absent could be related to a barrier-overcoming or tunneling for shear-thinning matter 
along the wavy-roughness (geometric valley and peak served as atomic potential surfaces) in 
annular tubes when the wavy-roughness is present. Once the geometry-tuned potentials (energy) 
overcome this barrier, then the tunneling (spontaneous transport) inside wavy-rough annular 
tubes occurs. 

Finally, We demonstrate in Fig. 1 that if we select the activation energy to be 4 x 1CT 19 J we 
can then observe a sudden drop of the resistance (frictional or shear stress) around 3 order of 
magnitude at T = 300. 5°K (Vh ~ 3.12 x 10~ 21 m 3 ). It means there is a rather low resistance 
below the temperature : T ~ 300°K for the material parameters selected. As r ~ (below 
T ~ 300 K), from |r| = V25Q/2 (rz 7^ 0), we can understand that there is no need for any 
external (gravity or electric field) forcing (5Q ~ 0) once the persistent current occurs. 
The possible reasoning for this frictionless transport of many condensed (composite) fermions 
can be stressed again as a brief summary. It could be due to the strong shearing driven by 
larger external (say, gravity or electric field) forcings along a confined wavy-rough tube. The 
shear-thinning (the viscosity diminishes with increasing shear rate) reduces the viscosity for the 
transport of this condensed many-fermion system significantly. One possible outcome for almost 
vanishing viscosity is the nearly frictionless transport. We shall investigate other relevant issues 
[18-20] in the future. 
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Fig. 1 Calculated (shear) stresses or resistance using an activation energy 4 x 10 19 J. 
There is a sharp decrease of shear stress around T ~ 300. 5° K. Below around 300 K 
(V h « 3.12 x 10~ 21 m 3 ), the transport of many composite (condensed) fermions is nearly 
frictionless. 



